Mass of nonrelativistic meson from leading twist distribution amplitudes. 
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In this paper distribution amplitudes of pseudoscalar and vector nonrelativistic mesons are con- 
sidered. Using equations of motion for the distribution amplitudes, it is derived relations which 
allow one to calculate the masses of nonrelativistic pseudoscalar and vector meson if the leading 
twist distribution amplitudes are known. These relations can be also rewritten as relations between 
the masses of nonrelativistic mesons and infinite series of QCD operators, what can be considered 
as an exact version of Gremm-Kapustin relation in NRQCD. 
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I. INTRODUCTION 

o : 

Study of hard exclusive processes can lead to a better understanding of mechanisms of hadrons production and 
Q , hardrons structure. The description of such processes is based on the factorization theorem [lj, Within this 
theorem the amplitude of hard exclusive process can be separated into two parts. The first part is partons production 
^1 | at very small distances, which can be treated within perturbative QCD. The second part is hardronization of the 
partons at larger distances. This part contains information about nonperturbative dynamic of strong interaction. For 
hard exclusive processes it can be parameterized by process independent distribution amplitudes (DA) , which can be 
r— i considered as hadrons' wave functions at light-like separation between the partons in the hadron. 

' Recently, the leading twist (twist- 2) DAs of nonrelativistic mesons have become the object of intensive study 
^ d, 0, i, H, 0, H, H, E3, 0, E3, ESI ■ In papers 0, [H[ it was derived relations between twist-2 and twist-3 DAs of light 
O : mesons. Unfortunately, these relations contain DAs of higher Fock states, for which one has very poor knowledge. 
For nonrelativistic mesons, there states are suppressed by higher powers of relative velocoty of quark- antiqaurk pair 
inside meson what leads to simplification of these relations. In papers [l?], Ell it was shown that if one ignores 
the contribution of higher Fock states, one can write the solution of equations of motion for twist-3 DAs as an infinite 
t-H ■ series in Gegenbauer polynomials. Actually, using equations of motion one can find exact formula for the twist-3 
DAs of nonrelativistic mesons through the leading twist DAs. In this paper we apply this exact formula and derive a 
relation between the masses of pseudoscalar and vector nonrelativistic mesons and leading twist DAs of these mesons. 

This paper is organized as follows. Next section is devoted to the derivation of the relation between mass and the 
leading twist DA of pseudoscalar nonrelativistic meson. In section three similar relation will be derived for the vector 
meson. 
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II. PSEUDOSCALAR NONRELATIVISTIC MESONS. 

In this section the relation which allows one find the mass of pseudoscalar nonrelativistic meson if the leading twist 
DA of this meson is known will be derived. To do this one needs the following DAs 

(P(p)\Q{z) lal5 [z-z]Q(-z)\0)^ = if PPa [ d^W^M), 



M JL_ f w.<, '■'>:> 



2M Q {n) 



{P[p)\Q{z) l5 [z,-z]Q{-z)\Q)^if P — / d£ e H*-VH£,M), (1) 



l 



where fi is the scale at which the axial and pseudoscalar DAs </?a(£, m)j <Pp(£> m) are defined, Mq(/i) is the running 
mass of the heavy quark Q, Mp is the mass of the pseudoscalar meson P, £ is the fraction of the whole momentum 
of the meson P carried by quark-antiquark pair, [z, —z] is a gluon string which makes the corresponding operators 
gauge invariant. 

Now few comments are in order. The meson P is assumed to be nonrelativistic. The width of DAs of nonrelativistic 
meson is of order of £ 2 ~ v 2 , where v is relative velocity of quark-antiquark pair in the meson P [5|. In the end point 
region |£| ~ 1 the motion of quark-antiquark pair is relativistic, so DAs in these regions must be strongly suppressed. 
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One can even assume that DAs of nonrelativistic meson are zero in the end point region. However, as was shown in 
paper || radiative corrections, which can be taken into account through the evolution equation for the DAs, lead to 
the nonzero effect of the DAs in the end point region. In this paper all DAs and all scale dependent quantities will 
be taken at the scale at which radiative corrections are not very important and below explicit scale dependence will 
be omitted. For instance, in paper [5| it was shown that for charmonia mesons the effect of the evolution of the DAs 
is negligible at scale \i ~ fio = 1.2 GeV. It should be noted that at this scale ~ /xo the DAs ¥?a(£, /•*)> V-p(£> A 4 ) are 
equal to each other at the leading order approximation in relative velocity [5[ . 

In paper [lij it was shown that applying equations of motion one can derive the following relations between the 
moments of the DAs </2a(£), </3p(£) 

(Dp = 5 n0 + ^1 ((C- 2 )p - r{e~ 2 ) A ) , (2) 



where r = 4Mq/M 2 , the moments (£™)a,p are defined as 

(Dap = / d£ £> A ,p(£). (3) 



i-i 



It should be noted that actually equations ^ contan also terms which describe the contribution of higher Fock state 
quark-antiquark-gluon. However, for the nonrelativistic mesons such states are suppressed by large power of small 
relative velocity v of quark-antiquark pair . For this reason in equations ^ and below such terms are disregarded. 
Equations ([2]) can be easily solved. The solution is given by the following expression 

MO = C + rJ^^dt, (4) 

where the constant C can be determined from the normalization condition for the function y>p(£) 

1 d£ ip P (0 = 1 = 2C + r £ d£ * + ? <p A (t). (5) 

Note that equation (TJJ unambiguously determines the twist-3 DA </3p(£) through the twist-2 DA (pa(Q- Actually this 
can always be done if one disregards the contribution of higher Fock states. One more example which demonstrates 
the last statement will be given below for vector nonrclativisitc meson. 

As was explained above the DA <pp(£) is zero when £ is sufficiently near to the end points £ — > ±1. In order to 
obtain this behaviour when £ — > ±1 for the function (fp(£) in equation (T?)) one must set C = 0, since the integral in 
the right hand side of equation (TJJ is zero for £ ~ ±1. If one puts C = to equation ([5]), one gets the formula 

Mp f 1 j, 1 + £ 



' ¥M(£). (6) 



4M§ j_ x ^(i-e) 

The running quark mass Mq here is taken at the scale at which radiative correction can be disregarded. Formula (J5]) 
is one of the main results of this paper. 

It is known that the moments (£")yi can be expressed through the QCD operators 

(P(p)\Qlu7 5 (iz"D a ) n Q\0)=ifpp 1/ (zp) n (C)A- (7) 

Here 

D = D-D, D = d -igB a (X a /2). (8) 

From this perspective, equation ^ is exact expansion of the mass of the pseudoscalar meson P in a infinite series of 
QCD operators (£") - (P{p)\Q lu j 5 {iz°D a ) n Q\Q) 

M 2 °° 

—£ J = Y,( n + !)(n = 1 + 3(£ 2 >a + 5(£ 4 )^ + 7(e) a + 0((e) a)- (9) 

Q n=0 



Recall that equation ^ is valid at the scale when radiative corrections are not very important. 
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Formula ([9]) has very simple interpretation. To understand what does it mean let us recall that the moments 
has the meaning of the average of the n-th power relative velocity of quark- antiquark pair in the direction of meson 
motion @ 1 . It is not difficult to understand that to get the average of the n-th power relative velocity (v n ), the 
(w™) must be multiplied by a factor of n+ 1, since to get average one should take the integral ~ J dcos^cos" 9. So, 
one has the series 1 + (v 2 ) + (v 4 ) + (v 6 ) + which can be symbolically written as 

M 2 1 

" =(P(p)\——\0). (10) 



4M, 



Evidently, this is quantum expression for the energy E — 2M/yl — v 2 of two relativistic particle with zero net 
momentum mov ing with the speed v. 

In papers [13, d, Eil it was pointed to the duality between light cone formalism and NRQCD description of hard 
exclusive processes. Although it was proved for some exclusive processes, today there is no strict proof of this duality. 
It would be interesting to know if it is possible to find formula in NRQCD dual to formula ©. To do this let us recall 
that in NRQCD one expands the mass of the pseudoscalar meson in operators 

,..nx_ 1 ™|X + (5)>|0) 



* ' M 2 (P(p)\ X + <p\0) ' 

where tp and x are bispinors which create antiquark and destroy quark respectively. To the leading order approximation 
one has 3(£ 2 ) = (v 2 ) |5|]. If one puts this equation to ((9]) one gets 

M 2 

_^ = l + ( ^ )+0((u 4 }) _ (12) 

which is well known Gremm-Kapustin relation [2(J. So, formulas ©, © can be considered as an exact version of 
Gremm-Kapustin relation. 

Equation ^ can be used to test models of DAs. As an example let us take the models of IS and 2S charmonia 
mesons proposed in papers 5, 7] at the central values of the parameters of these models. For the running mass of the 

c-quark M c the value M c = M^(fi = AlJ^) = 1.23 will be taken [21]. Thus one gets M Vc = 2.81 GeV, M„/ = 3.53 
GeV, which is rather good estimation. 

III. VECTOR NONRELATIVISTIC MESONS. 

In this section the formula similar to formula © but for the vector nonrelativistic meson will be derived. To do 
this one needs the following DAs 

(V(p, e)\Q a (z)Q p (-z)\0) = J 1 d^)i | fe - Vx(0 +pS^(0 + 



h 1 , v h 2M Q 



fy My ^ > + 8 I 1 - fv M?) WrW*' MO)^, (13) 

where the constants /r, fv are defined as 

(V(p,e)\Q la Q\0) = f v M v e a , (V(p, e)\Qa aP Q\0) = ifr(p a ep - p e a ). (14) 

The DAs <pl(Q, Pt(Q are twist- 2 DAs, the DAs <£>_l(£), <Pa{£,) are twist-3 DAs. The equations of motion for DAs (fH?)) 
give [II] 

(n + 1)(D± = (Di + fcil (1 - S)(C- 2 )a, (15) 
\(n + 2)(l-S)(C)A = (C)±-S(C)T, 



1 It is assumed that the meson is moving along z-direction. 
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where (C 1 )l,t,±,a are the moments of the DAs (Pl{.£), <£t(£)> <P-L-(£)i Va(€), <Km) = 2 /t/ fyM Q /My. 
Equations (TiTi|) can be exactly solved [la]. For instance, the solution for the DA cp^ is 

^(0 = c + de + ^ (c 1 + / TT7^ + o- - o /* r^*) ■ ( 16 ) 

where — Vi(0 + ^£Vt(0- The function (^(C) must be even, normalized and it must be suppressed in the end 
point region. These requirements are fulfilled if Co = C\ = and 8 is equal to 

/t2Mq_ J^iijr^O 

This formula is the analog of equation ([6]) but for vector meson V. It is more complicated than equation ([6]) since 
vector meson has two independent constants fx and fy . Note that if one expands now right hand side of the last 
equation and puts Gremm-Kapustin relation (Tl~12|) one will get NRQCD expression for the ratio fr/fv [22[ |23| 

f = l-^ (18) 
fv 6 
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